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ABSTRACT 

We present results from a statistical study of clouds in two-dimensional numerical 
simulations of the interstellar medium. Clouds are defined as connected sets of pixels 
above an arbitrary density threshold. Similarly to real interstellar clouds, the clouds 
in the simulations exhibit a differential mass spectrum dN(M)/dM ~ M~ 1A4±A and 
a velocity dispersion-size relation Av ~ _R°- 41± - 08 ; although the latter with correlation 
coefficients of only ~ 0.6. However, the clouds do not exhibit a clear density-size relation. 
At a given mean density, clouds span a range of sizes from the smallest resolved scales 
up to a maximum given by a Larson-type relation -R max ~ p a , with a = —0.81 ± .15, 
although numerical effects cannot be ruled out as responsible for the latter correlation. 
The continuum of sizes at a given mean density implies a range of column densities of 
up to two orders of magnitude. 

This result supports the suggestion by Scalo that the observed relation may be an 
artifact of the observational limitations. In this case, the non-existence of a density-size 
relation would suggest that the origin of the Av-R relation is not virial equilibrium of 
clouds that follow a p oc R^ 1 law. Instead, the Av-R relation can be interpreted as 
a direct consequence of a k~ 2 turbulent spectrum, which is characteristic of a field of 
shocks, verified in the simulations. However, we also discuss the possibility that the 
clouds are in balance between self-gravity and turbulence, but with a scatter of at least 
a factor of 10 in the velocity dispersion-size relation, and a scatter of over a factor of 
100 in the density-size relation, according to the equilibrium relation Av ~ (NR) 1 ^ 2 . In 
this case, the scatter in column density would be larger than the dynamic range of the 
simulations (and most observational studies). 

Finally, we compare these results with recent observational data. We propose a 
simple model suggesting that recent results in which nearly constant column densities 
are derived for dark IRAS clouds may be an artifact of a temperature gradient within 
the clouds induced by external radiative heating. As a consequence, we emphasize that 
IRAS surface brightness maps are not appropriate for measuring column densities. 



Subject headings: hydrodynamics — ISM: clouds — ISM: structure — turbulence 
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1. Introduction 

Interstellar clouds appear to follow a set of scal- 
ing relations first noticed by Larson (1981), and then 
apparently confirmed (with slight modifications) by a 
number of other workers. These "Larson's relations" 
have the form 



p~R a 
Av ~ R p , 



(la) 
(lb) 



where R is the cloud size, p is the gas density, Av is 
the velocity dispersion derived from the line widths, 
and a and (3 are the constant scaling exponents. Ad- 
ditionally, the clouds are found to exhibit a mass dis- 
tribution of the form 



dN(M) 
dM 



M n . 



(2) 



The most commonly quoted values of the exponents 
area 1.15±.15, (3 ~ 0.4±.l, and n 1.55± .15 



(Larson 1981 




Lorrelles et al. 1983; 


Dame et al. 198C 


Falgaronc & Pcrault 1987; Myers & Goodman 1986 



Falgarone, Puget fe Perault 1992, hereafter FPP 



Fuller fc Myers 1992; Miesch fc Bally 1994; Wood, 



Myers fc Daugherty 1994, hereafter WMD; Caselli fc 



Myers 1995| see also the reviews by Scalo (1985, 1987) 
and Blitz (1991)). However, significantly discrepant 



values have also been reported (e.g., Carr 1987; Lorcn 



1989), and the validity of these scaling relations is 
currently the subject of strong controversy within the 
community. 

The above "standard" values of the exponents for 



of the virial theorem (e.g., 


Larson 1981 


; Myers & 


Goodman 1988; 


Caselli & Myers 1995). For a = -1 



(which coincidentally implies constant column den- 
sity), a value (3 — 0.5 implies virial balance between 
self-gravity and the internal velocity dispersion. How- 
ever, note that for an arbitrary value of the density 
scaling exponent a, a correspon ding virial balance 
value of (3 can always be found ( |Vazqucz-Scmadeni 
'k Gazol 1995 ). Thus, the density-size relation (|la|) 
remains unexplained. In fact, it has been proposed 
by Scalo (1990) that this relation may be a mere ar- 
tifact of the dynamic range limitations of the obser- 
vations, and does not reflect a real property of inter- 
stellar clouds. In particular, in the case of molecu- 
lar line data, the observations are restricted to col- 
umn densities large enough that the tracer molecule 
is shielded against photo-dissociating radiation. On 



the other hand, while the proportionality between 
line integrated CO intensity and mass surface density 
has been reliably established for extragalactic obser- 
vations (Dickman, Snell, & Schloerb 1986), this rela- 



tionship is only valid for scales at which calibrations 
have been possible, i.e., scales larger than a few pc. 
Furthermore, for clumps within molecular clouds, the 
structures identified in CO often do not correspond 
to those identified with higher-density tracers (e.g., 
Massi & Lizano 1994; J. Scalo, private communica- 



tion). 

In this paper we present the mass spectrum of 
clouds and search for Larson-type correlations in 
three two-dimensional numerical simulations of tur- 
bulence in the interstellar medium (ISM), one from 
Passot, Vazquez-Semadeni & Pouquet (1995, here- 
after Paper I), and the other two being variants of 
the former with respectively larger density contrasts 
and larger resolution. We have developed a cloud- 
identifying algorithm which allows us to measure the 
average density, velocity dispersion and total mass 
within well-defined (though arbitrary) cloud bound- 
aries in the density fields of the simulations. The 
cloud sample thus obtained has the advantage over 
actual observations that no tracer limitations exist in 
"detecting" the clouds, although, on the other hand, 
a form of saturation is present due to the relatively 
small density dynamic range and other numerical lim- 
itations of the simulations. 

In § 2 we briefly describe the simulations and the 
cloud-identifying algorithm, and in § 3 we present the 
resulting statistical cloud properties. In § 4 we ana- 
lyze the limitations of the simulations, discuss the im- 
plications of the absence of a density-size relation, and 
give a preliminary discussion of the physical mecha- 
nisms behind the velocity dispersion-size relation. We 
also compare the results with recent corresponding 
observations, in particular those of WMD. Finally, § 
5 summarizes the results. 

2. Numerical Method 

Paper I presented magneto-hydrodynamical simu- 
lations of the ISM incorporating model terms for cool- 
ing, diffuse heating and local heating from star for- 
mation. The simulations solve the equations for the 
evolution of the density, velocity, internal energy and 
magnetic fields in the presence of self-gravity, namely 



^ + V • (pa) = pV 2 P , 



(3) 



2 
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at p 

(5) 

— = V x (u x B) - j/ 8 V 8 B, (6) 



p-1. 



(J) 



We refer the reader to Vazquez-Semadeni, Passot & 
Pouquet (1995a) and Paper I for full details on the 
simulations. Here we just point out that the simu- 
lations from which the data are extracted represent 
a square section of the ISM along the Galactic plane 
of size 1 kpc on a side, with a resolution of 512 grid 
points per dimension. Also, it is important to note 
that all the evolution equations contain dissipative or 
diffusive terms which are necessary since the numer- 
ical technique used to solve the equations (spectral 
method) does not produce numerical viscosity, so the 
dissipation must be included explicitly. The momen- 
tum and magnetic field equations contain "hypervis- 
cosity" terms of the form V 8 , which confine the vis- 
cous effects to the very smallest scales in the simula- 
tions ( [Babiano et al. 1987| [Mc Williams 1984[ ). How- 
ever, although globally dissipative, hyperviscosity is 



not everywhere positive definite (Passot & Pouquet 
19^8|), so standard Laplacian terms are used in the 
continuity and internal energy equations. 

Throughout the paper, densities are expressed in 
units of 1 cm -3 and velocities in units of 11.7 km s _1 , 



the units used in Paper I. In particular, in the present 
pap er we will use data from the run labeled Run 28 



pam 
in E 



aper I , However , the star- formation scheme used 



in the simulations of Paper I assumes that a star is 
formed wherever the local density exceeds a critical 
value p c . (A "star" in the simulations is a point source 
of heat.) This naturally imposes an upper limit on the 
densities reached by the model, since the stellar heat- 
ing increases the local pressure and causes the gas 
to expand. Thus, densities above p c are very rarely 
reached. In Run 28, p c = 30, limiting the density con- 
trast Pxaax/ Pmin to values ~ 1000. In order to obtain 
a somewhat larger dynamic range, we have performed 
an additional run, called Run 28bis, which is identical 
to Run 28 up to t — 6.5 x 10 7 yr, but afterwards has 
the star formation turned off. This run ends up col- 
lapsing gravitationally at t ~ 8.8 x 10 7 yr due to the 



lack of support from stellar energy injection, but in- 
termediate times provide a good framework for study, 
exhibiting density maxima p> 100 cm~ 3 , and density 
contrasts J> 5000. Finally, a run similar to Run 28bis 
but at a larger resolution (800 x 800, refered to as 
Run 28.800) was also performed to discuss the effects 
of dissipation. As an illustration, fig. [I] shows a con- 
tour plot of the density field of Run 28 at t = 6.6 x 10 7 
yr, with the velocity field represented by the arrows. 
Note that, although the size scales represented by the 
runs are larger than those studied by most papers con- 
cerned with cloud statistics, we expect the results to 
be applicable, since all terms in the equations solved 
in the simulations are scale-free, except for the dissi- 
pative ones. 

In order to investigate the statistical properties 
of the clouds in the simulations, we have developed 
an automated algorithm which identifies and labels 
clouds. We define clouds as connected sets of pixels 
with densities above a user-defined density threshold 
p t . The types of clouds that are identified in this way 
are thus clearly dependent on the value chosen for p t , 
a somewhat similar situation to performing observa- 
tions using various density tracers. As an example, 
figs. ||a and |^b show the clouds obtained by respec- 
tively setting p t = 4 and p t — 16, at time t = 6.6 x 10 7 
yr into the evolution of Run 28. 

The extremely irregular shapes of most clouds are 
noteworthy, and in fact pose a problem in defining 
the "size" R of a given cloud. For simplicity, here 
we take the size of a cloud as the square root of the 
number of pixels it contains. This definition may be 
somewhat unrealistic if clouds are fractal in such a 
way that their perimeters are not proportional to the 
square root of their areas (e.g., |5calo 1990 ; Falgaronc 



Phillips & Walker 1991), but for simplicity and sim- 
ilarity to observational work we adopt it here. Also, 
for convenience, sizes are expressed in pixels in most 
of the figures below (1 pixel ~ 2 pc). 

For each cloud found by the algorithm, it is then 
a trivial matter to compute the average density, the 
mass, and the turbulent velocity dispersion, calcu- 
lated as At; = ((u — (u)) 2 } 1 / 2 , where u is the local 
value of the velocity and (u) is the average over the 
cloud area. Note that this quantity is not density- 
weighted. 
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3. Statistical Cloud Properties 

We have obtained relatively large samples of clouds 
by considering values of p t =3, 4, 6, 8, 11, 16, 23, and 
29 in Run 28, and p t =3, 4, 6, 8, 11, 16, 23, 32, 45, 64, 
and 90 in Run 28bis. Clouds obtained with each value 
of pt are indicated with a particular symbol type in 
the plots discussed below. In what follows, the data 
from Run 28 will be considered at t = 6.6 x 10 7 yr, 
and Run 28bis at t = 7.15 X 10 7 yr, unless otherwise 
noted. The cloud samples for those cases contain 158 
and 145 clouds, respectively. 

Figures ||a and |^b show the logarithmic mass dis- 
tributions of the clouds for Run28 and Run 28bis, 
respectively. The masses are expressed in the nondi- 
mensional units of the simulations, in which the total 
mass contained is An 2 , corresponding to 0.36 x 10 7 Mq 
in real units. The turnover at low masses may be 
attributed to incompleteness due to the finite resolu- 
tion. The high-mass sides of the distributions, how- 
ever, exhibit least-squares slopes (fitted for logM > 
0.8 and shown as the straight lines) which imply in- 
dices n — —1.51 and —1.43, respectively (c.f. eq. [||). 
These values are within the range of values reported 



in various observational results (e.g., WMD; see also 
the reviews by |Scalo 1985| , |Blitz 1991| , |Mundy 199l| , 
and references therein). 

Figures ||a and |]b show the correlation between 
velocity dispersion and size for all clouds with sizes 
J> 2 pixels (1-pixel clouds are excluded as they have 
no velocity dispersion). Least squares fits to the data 
give slopes (shown as the solid lines) (3 ~ 0.37 and 
f3 ~ 0.39 respectively for Run 28 and Run 28bis, with 
moderately tight correlation coefficients ~ 0.6. The 
large scatter of about one order of magnitude in the 
correlations should be emphasized. This is compara- 
blcto the scatter found for molecular- line data (e.g., 
FTP). The derived values of j3 <~ 0.4 are slightly lower 



than the most commonly accepted value of 0.5 (dot- 
ted lines), but are remarkably close to d eterminations 
that i nclude heterogeneou s data sets ( Larson 198l| ; 



FTP; Fuller fc Myers 1992) 



Figures ^|a and[^b show the average density (p) vs. 
size for the snapshots of the two runs. Several points 
should be remarked. First, note that, in general, the 
average density of the clouds is quite similar to the 
value of pt used to define the cloud, and in fact tends 
to increase with size at each p t . This is not surprising, 
due to the presence of dense "cores" inside the largest 
clouds at each p t , which tend to increase their average 



density. Second, no clear correlation can be seen in ei- 
ther figure. Instead, at a given p t , clouds down to the 
smallest possible size are seen. These are small clouds 
with low densities and, therefore, low column densi- 
ties^] N = pR. Third, the size of the largest cloud at 
each pt is smaller for larger p t . In particular, in both 
figures the set of largest clouds at each p t seems to lie 
near ap~ R^ 1 law, similar to the standard exponent 
in Larson's relation ((Taj). 

In order to test the robustness of the above results, 
we performed the same analysis at various other times 
in both runs, namely at t =7.8, 9.1, 10.4, 11.7 and 
13.0xl0 7 yr in Run 28, and t =7.8 and 8.6xl0 7 yr in 
Run 28bis. The average value of the exponent /3 of 
the velocity dispersion-size relation is ((3) = 0.41 ±.08, 
with typical correlation coefficients ~ 0.6. For the 
mass spectrum, an average exponent (n) = — 1.44 ± .1 
is found. The errors are the standard deviations of the 
set of values found for all times. These results con- 
firm the fact that the simulations show correlations 
consistent with the observations in both cases. 

Regarding the density, in all cases small clouds 
with low densities exist, the plots (not shown) be- 
ing qualitatively similar to figs. ||a and ^Jd, and the 
full ensemble of clouds not exhibiting any correlation 
with size. The largest clouds at each p t , on the other 
hand, continue to exhibit a near power-law relation 
with size, with average exponent (a) = —0.81 ± .16. 
This is smaller than the slopes found in figs. ^|a and 
||b, which coincidentally seem to have some of the 
steepest slopes in the distribution. This is illustrated 
in fig. ^, which shows density vs. size for the largest 
cloud at each value of pt at all times considered above 
for the two runs. In this figure, clouds obtained with 
a given value of p t at any one time in cither run are 
shown with the same symbol and joined by a dotted 
line. The resulting curves have been displaced by in- 
crements of 0.2 in log(p) for clarity For reference, 
the solid line shows a (p) oc R^ 1 power law. The 
three uppermost curves in the figure correspond to 
Run 28bis and, because of the larger density dynamic 
range of this run, exhibit power-law behavior over a 
larger range of scales, while those for Run 28 saturate 



1 Note that the column density defined here refers to a cut 
through the clouds on the plane of the simulations, and has 
units of cm -2 . To obtain a column density with units cm -1 , 
appropriate to our two-dimensional problem, a multiplication 
by the unit length along the third (2) dimension, perpendicu- 
lar to the plane of the simulation, should be performed. For 
simplicity, we omit this constant factor throughout the paper. 
A similar situation applies to the computation of masses. 
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at p ~ 30. For this reason, only clouds with logi? > 1 
for Run 28bis, and with logi? > 1.5 for Run 28 were 
considered in computing (a). The implications of 



these results are discussed in 4.2.1 For convenience, 



in what follows we drop the brackets when referring 
to the average density of clouds. 

Finally, it should be pointed out that, although 
clouds with low column densities exist in the simu- 
lations, most of the mass still resides in the largest 
clouds, since the distribution of cloud sizes at a given 
mean density appears to be roughly uniform. How- 
ever, this is possibly an effect of the absence of super- 
novae in the simulations, since the "expanding HII 
regions" included do not have enough power to dis- 
rupt the largest gravitationally bound complexes. 

4. Discussion 

4.1. Applicability and limitations of the re- 
sults 

The results presented in § || have important im- 
plications, provided that the simulations are indeed 
representative of ISM dynamics. That this is likely 
to be the case, in spite of their two-dimensionality, is 
suggested by the fact that the simulated ISM repro- 
duces both the velocity dispersion-size relation and 
the mass spectrum of the clouds (c.f. § ||), as well as 
other physical properties of the ISM, such as the mean 
density of giant cloud complexes flVazqucz-Scmadeni 
it al. 1995a), the cloud and intercloud magnetic field 



strengths (Paper I), the rate of formation of massive 



stars ( Vazquez-Semadeni, Passot & Pouquet 1995b), 
etc. However, one important possible criticism due 
to the two-dimensionality, is that a p ~ Rr 1 scal- 
ing relation in three dimensions might translate into 
a p ~ R° relation in two dimensions, just because of 
the elimination of one dimension. This could be at the 
origin of the near constancy of (p) observed at every 
value of pt ■ Closer examination shows that this argu- 
ment is invalid. Even though the clouds seem to have 
nearly constant densities at each value of pt , this only 
reflects the fact that, due to our cloud-identifying al- 
gorithm, small, dense clouds are not seen at small p t , 
since they are "hidden" within the larger, low density 
ones. But obviously the density is not constant in 
the simulations, invalidating the possibility of an R° 
dependence of the density. 

Another important source of concern is that a sig- 
nificant fraction of the clouds in the samples have 
sizes of only a few pixels, and their properties are 



thus likely to be affected by viscosity and diffusion. 
Thus, it is important to quantify the extent to which 
the results of the previous section might be influenced 
by these terms. In particular, the question arises as 
to whether the existence of low-column density, small 
clouds might be an artifact of the dissipative terms. 

Concerning viscosity, owing to the hyperviscosity 
scheme with a V 8 operator, its effects on the veloc- 
ity field are confined to wavenumbers k in the range 
fcmax/2 <; k < fc max with fc max = 255 for Runs 28 
and 28bis, as can be seen from the spectrum in fig. 
HO (see §4.2.2). The same applies to the dissipative 



term in the magnetic equation (see fig. 5 of Paper I). 
Naively, one would then expect the effects of viscosity 
to be present at scales up to twice the smallest scale 
of the simulation, i.e., from one to two pixels. Actu- 
ally, the correspondence between scale ranges in real 
and Fourier spaces is not as sharp, and one can expect 
"leakage" up to possibly 4 pixels. Visual inspection 
of the velocity field confirms that shocks are spread 
over typically 4 pixels. Nevertheless, note that vis- 
cosity is not effective if the velocity gradients are not 
large, and thus the 4-pixel estimate is an upper limit 
to the sizes affected by viscosity. It is worth pointing 
out that in Paper I the range of influence of viscos- 
ity was estimated at ~ 5 pixels. However, this was 
an over-conservative estimate not based on a detailed 
analysis of shock widths, and the 4-pixel figure given 
here should be considered as a slightly more precise 
estimate based on the above considerations. 

In order to correct for viscosity effects, clouds with 
sizes up to 4 pixels in size are excluded from the "cor- 
rected" plots shown in the Appendix. Note that, since 
clouds have in general elongated shapes and sizes are 
computed as the square root of the number of pixels, 
the possibility exists that clouds with computed sizes 
larger than 4 pixels will still be 4 pixels or less across 
one particular direction. However, we believe this ef- 
fect may be roughly compensated by the fact that the 
4-pixel estimate is an upper limit, and thus we retain 
all clouds with sizes larger than this. 

Of greater concern are the possible effects of mass 
diffusion, since the standard Laplacian diffusive oper- 
ator used in the continuity equation causes diffusion 
to be important over a larger range of scales than 
the hyperviscosity. Indeed, the characteristic diffu- 
sion time Tdif can be shorter than the turbulent cross- 
ing time tnl for clouds smaller than about 16 pixels 
(see Appendix), and diffusion may dominate over tur- 
bulent advection for those clouds. In particular, it is 
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possible that the small, low-density clouds reported 
in the previous section might be a numerical artifact 
of the diffusion, which tends to reduce density peaks 
and spread the clouds out, or, conversely, to prevent 
clouds from reaching higher densities and smaller sizes 
than they do in the simulations. 

Note, however, that the effect of mass diffusion is 
exclusively to damp density gradients originated by 
the turbulent velocity field, and so diffusion is in- 
capable of forming clouds by itself; instead, it only 
modifies the properties of clouds formed by turbu- 
lence or other processes (gravity or the various insta- 
bilities discussed in Paper I). We can thus obtain a 
crude estimate of the size and maximum density a 
cloud would have in the absence of diffusion by inte- 
grating the diffusion equation backwards in time for 
initial conditions corresponding to the clouds in the 
simulation, over the length of a nonlinear time tnl- 
This correction overestimates the effects of diffusion, 
since it neglects the advection term entirely, and the 
nonlinear time is computed using the turbulent ve- 
locity associated to the size of the cloud, as given by 
the turbulent spectrum in fig. [lo[ However, it is actu- 
ally turbulent scales larger than a cloud's size which 



are most likely to form it (Elmegreen 199S ; Vazquez- 
Scrjiadeni, Passot fc Pouquct 1995b| ). In particular 



in the simulations at least, clouds often form at the 
interfaces of expanding shells from previous star for- 
mation events, which may have velocities of several 
kilometers per second, rather than the significantly 
smaller velocities (<; 1 km s _1 ) indicated by the spec- 
trum, which is a globally averaged quantity. Nev- 
ertheless, for robustness, we will use the worst-case 
correction. The details of this calculation are given 
in the Appendix. There it is shown that the cen- 
tral density of a cloud varies by factors of 3-5 in the 
worst cases under the influence of diffusion. Using 
this "correction" factor and assuming the clouds move 
along lines of constant mass in the density-size plot, 
one can produce a "corrected" such diagram, shown 
in fig. |l6| in the Appendix for Run 28bis. There it 
can be seen that, although clouds are indeed brought 
slightly closer to a correlation, at the lowest average 
densities cloud sizes still vary by factors of about 100, 
maintaining the conclusions from 

As a further test, we have also produced a pre- 
liminary higher-resolution run, labeled Run 28.800, 
similar to Run 28bis, but at a resolution of 800 x 800 
pixels, the largest that we can perform in the CRAY 
YMP of DGSCA, UNAM. The density field for this 



run is shown in fig. ^. In order to produce this run, 
the data from Run 28 at t = 6.5 x 10 7 yr were inter- 
polated to produce initial conditions for the 800 x 800 
simulation, and then evolved for another 0.65 x 10 7 yr, 
enough to develop the additional small scale structure 
corresponding to the larger resolution. Incidentally, 
it is worth noting that, even though Run 28.800 has a 
resolution only ^60% larger than Runs 28 and 28bis, 
the computational effort it requires is roughly 5 times 
larger in CPU run time. Fig. || shows the resulting 
density-size raw plot for this run, including cloud sizes 
down to 1 pixel. In the Appendix, fig. |l7| shows the 
corresponding plot incorporating the corrections de- 
scribed above, namely elimination of clouds with sizes 
< 4 pixels, and the correction for diffusion. In order to 
maximize the available dynamic range, in the figures 
in the Appendix we have used values of the density 
threshold as low as p t = 0.5. For Run 28.800, it is 
seen that the clouds with the lowest average densities 
span a range of roughly a factor of 200 in sizes after 
the corrections, while the raw data exhibit a range of 
a factor ~ 500. 

In summary, the discussions from this section sug- 
gest that the existence of small clouds with low den- 
sities described in §|jj| is a true consecuence of the dy- 
namics and not an artifact of the dissipative terms 
used to stabilize the equations. 

4.2. Implications of the results 

4-- 2.1. Density, size and equilibrium 

With the above considerations in mind, we then 
have a number of direct implications. First, as already 
stated above, the absence of a density-size relation 
implies non-constant column densities for the clouds. 
Specifically, cloud column densities vary by over two 
orders of magnitude in the simulations. Additionally, 
the result that clouds with sizes down to the smallest 
scales exist at all values of the mean density implies 
that the observed density-size relation may indeed be 
a product of the limited dynamic range of typical 
surveys. Under these conditions, the derived column 
densities could appear to show variation over a mod- 
est range only, while the observed sizes might span a 
range larger than three orders of magnitude, thus cre- 
ating an apparent correlation. Observational effects 
as a possible origin of the density-size relation were 
first suggested by Larson (1981) himself, and later 
discussed in more detail by Scalo (1990). Another 
reason for the appearance of a spurious density-size 
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relation may be a selection effect introduced by the 
tendency of observational work to focus primarily on 
global intensity maxima of the maps, therefore pos- 
sibly missing weaker, local maxima, a bias that only 
recently has started to be avoided (e.g., FPP). 

The Larson-type relation defined by the largest 
clouds, with a scaling exponent a ~ —0.8 ± .16 is 
particularly interesting. The immediate question that 
arises is whether this relation is physical, or is induced 
by numerical constraints. On the physical side, a first 
consideration is that, if the clouds are hierarchically 
nested (smaller, denser clouds are part of larger, less- 
dense ones), then mass conservation implies — a < 3 
in three dimensions. In our two-dimensional case, this 
limit becomes —a < 2. This is a physical limit which 
in the simulations is pushed closer to the observed 
relation because of the two-dimensionality. 

Another physical issue is whether the large clouds 
which follow a density-size scaling law are virialized, 
so that the standard scenario in which Larson's re- 
lations hold for virialized clouds would apply to the 
largest clouds. However, examining cloud virializa- 
tion in our simulations turns out not to be a straight- 
forward task(§4.2.2). Here, we just check whether the 
standard Av-R relation is also satisfied in the large 
clouds. To this end, some clouds in figs. ^ and || have 
been labeled with numbers so that they can be identi- 
fied from one figure to another. Interestingly, for Run 
28 at t = 6.6 x 10 7 yr (case a in the figures), the largest 
clouds, which are very close to the constant-column 
density line p ~ R , shown as the straight line in 
fig. 5 a, are also very close to the line Av ~ i? 1 / 2 in 
fig. fa, suggesting balance between gravity and tur- 
bulence for this set of clouds. However, this is not 
the case for the clouds in Run 28bis at t = 7.15 x 10 7 
yr (case b in the figures). In this case, although the 
clouds again show a slope very close to —1 (fig. ||b), it 
can be seen from fig. ^b that they all have compara- 
ble velocity dispersions. We conclude that even clouds 
with a p oc R^ 1 density dependence are not necessar- 
ily in equilibrium between turbulence and self-gravity. 
Possibly, magnetic support is more important for the 
latter set of clouds, as in the results of Myers & Good- 
man (1988). 

An interesting question is whether the various 
clouds along the large column density "envelope" of 
the distribution in the density-size relation are essen- 
tially the same cloud seen at various different thresh- 
olds p t , or else they are truly different clouds. In fact, 
the answer is that they are neither. This is exem- 



plified in fig. [| which shows a few selected branches 
of the cloud hierarchy for Run 28bis. One branch 
off the largest cloud includes all the clouds along 
the envelope, but both the largest and second largest 
clouds are seen to also have other branches to daugh- 
ter clouds off the envelope (dotted lines) . The same is 
true of clouds lying immediately below the envelope 
(dashed lines), which seem to define a second envelope 
of similar slope (see fig. ||) . 

On the numerical side, the mass diffusion term 
in the continuity equation may tend to reduce the 
column density of clouds defined through the thresh- 
old density-criterion we use here, since the diffusion 
widens and smooths clouds, whose outer parts may 
then be left out of the domain defined by p t . This 
effect, plus the plain limitations imposed by the res- 
olution, clearly prevent the formation of very small 
clouds with very large column densities, causing shal- 
lower slopes of the high-column density envelope. 
In fact, while the average slope we obtain implies 
smaller column densities for smaller clouds, it has 
been pointed out by Scalo (1985, sect. III. A) that it 
is obvious from inspection of Lynds' dark cloud cat- 
alog that smaller clouds are darker. Indeed, in our 
case, the correction discussed in 4.2.1 and in the Ap- 
pendix tends to bring this envelope towards steeper 
slopes. In summary, the specific slopes defined by the 
largest clouds in the logp-logi? plot cannot be un- 
ambiguously attributed to real physical effects. High- 
resolution, three-dimensional simulations are needed 
to resolve this issue. Unfortunately, the largest simu- 
lations of supersonic compressible turbulence known 
to us (e.g., Porter, Pouquet, & Woodward 1994, us- 
ing 512 3 grid points), are purely hydrodynamic and 
do not contain many essential ingredients of ISM dy- 
namics, such as the magnetic field, self-gravity, and 
stellar (i.e., small-scale, compressible) forcing. Thus, 
the necessary calculations are still a few years in the 
future. 

In any case, regardless of what the specific slope 
of the envelope turns out to be upon removal of nu- 
merical effects, our results suggest that the notion of 
a density-size scaling relation should probably be re- 
placed by that of an "allowed" region in p-R space. 
Whether the high-column density boundary is truly a 
power law, the value of its corresponding index, and 
the physical mechanisms responsible for it, are issues 
that remain to be determined by high resolution 3D 
simulations. 
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4-2.2. Velocity dispersion- size relation and turbu- 
lence 



spectral index m — — 2 (e.g., Bonazzola et al. 1987 



Vazquez-Semadeni & Gazol 1995; Padoan 1995; Fleck 



Since the density-size relation is not verified for the 

clouds in the simulations, yet the dispersion-size rela- 
tion 
rel 



is, the standard argument explaining the Av R 



tion, namely virial equilibrium in clouds satisfy- 
ing p oc cannot be invoked. This implication is 

independent of the dimensionality of the simulations, 
as it only relies on the non-existence of a density-size 
relation, and not on particular values of the scaling 
exponents. That is, the two-dimensionality would, at 
most, change the exponent in the virial equilibrium 
relation between Av and p, but not destroy the corre- 
lation altogether. Since in the simulations no density- 
size relation exists, a unique virial equilibrium rela- 
tion between Aw and p does not exist either (although 
see § 4.2.3| ). Note, however, that a Jeans- type analysis 



incorporating a "turbulent pressure" ( Chandrasekhar 
19SHJ 
195 



Scalo & Struck-Marcell 1984 




Bonazzola et al. 


Elmcgreen 1991 


; Vazquez-Semadeni & Gazol 

/A..\2r7. _• 



Av ~ Rp 



1/2 



(8) 



for clouds of size R equal to their Jeans length. Since 
the latter is independent of dimensionality, relation 
(||) holds also independently of dimensionality. Thus, 
the same scaling laws as in 3D are expected for two- 
dimensional clouds in balance between self-gravity 
and turbulence. This is consistent with a crude esti- 
mate of virial balance in which one equates the grav- 
itational energy W to the turbulent kinetic energy 
K. For simple cloud geometries and uniform densi- 
ties and turbulent velocity dispersions, in 3D one ob- 
tains GM 2 /R ~ M (Av) 2 . Taking M cx pR 3 gives the 
usual result (Av) 2 oc pR 2 . In 2D, on the other hand, 
M oc pR 2 , but the gravitational energy becomes 
W ~ GM 2 , thus preserving the result (Av) 2 oc pR 2 . 
However, such a simplified treatment may not be ap- 
plicable, since in 2D logarithmic corrections appear, 
and also the gravitational term in the virial theorem 
may differ from the gravitational energy. A detailed 



analysis of this problem is in progress ( Ballesteros- 
Paiedes & Vazquez-Semadeni 1996). Thus, the in- 



stability analysis is provisionally a more reliable indi- 
cator of the equilibrium relation expected in 2D. 

In the absence of a density-size relation, a plau- 
sible origin of the velocity dispersion-size relation is 
the statistical properties of the turbulence itself. In- 
deed, an index j3 = 1/2 is expected for turbulence 
characterized by an energy spectrum E(k) oc k m with 



For comparison, the straight 
2 power law. It is clearly seen 



1996), where k is the wavenumber associated with 
scale R — 2n/k. Such a spectral index is the sig- 
nature of a field of random shocks (see Passot, Pou-| 
quct fc Woodward 1988j and references therein) . Fig- 
ure |10| shows the spectra for the incompressible (solid 
line) and irrotational (dashed line) parts of the ve- 
locity fields of fig. [l]. 
dotted line shows a k 
that the spectrum of the incompressible component 
is remarkably well described by this slope. The irro- 
tational, or compressible, component exhibits some- 
what stronger fluctuations (most likely due to the fact 
that the "stars" in the simulations inject purely com- 
pressible energy), but is still very close on the average. 
The 20% discrepancy with the index /? ~ 0.4 found 
in the simulations may be due to the fact that in the 
simulations there exists an upper bound to the turbu- 
lent velocity dispersion (of order a few km s _1 ) that 
can be imparted to the medium by the stellar heating, 
since in the model they only heat the gas to ~ 10 4 K. 
This introduces a "truncation" on the Av-R relation, 
which flattens the resulting slope, as can be seen in 
figs. |a and |b. 

4-2.3. An alternative interpretation 

At this point, one important alternative must be 
pointed out. The large scatter in the velocity dispersion- 
size relation would be consistent with the clouds not 
having constant column densities, even if they were in 
equilibrium between turbulence and self-gravity. In- 
deed, in equilibrium, the scatter of about one order 
of magnitude in the velocity dispersion would imply 
a scatter of roughly two orders of magnitude in the 
column density, as can be seen from the equilibrium 
relation 

Av ~ (A^i?) 1/2 , (9) 

which is equivalent to relation (||). Thus, our results 
can also be interpreted in the sense that all clouds 
tend to be virialized, although with a scatter of at 
least two orders of magnitude in the column den- 
sity. Preliminary results on the energy budget of a 
smaller cloud sample in Ru n 28 ( Ballesteros-Paredcs| 
fc Vazquez-Semadeni 1995| ) suggest that the sum of 
the kinetic and magnetic energies is within one order 
of magnitude of the gravitational energy, although 
surface terms were not considered there. Also, the 
total gravitational and the turbulent kinetic energies 
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per unit mass in the simulations arc almost in equipar- 
tition, as shown in fig. |ll| 

In this alternative interpretation, a density-size 
correlation may still be present, but missed by the 
simulations because the scatter is larger than the col- 
umn density dynamic range in the simulations (and 
in most observational studies). This scenario would 
leave the origin of the putative density-size relation 
unexplained, although it may still be possible that 
the Av-R relation is originated by the turbulent en- 
ergy spectrum, and that the p-R relation is the con- 
sequence of virial equilibrium. Although this sce- 
nario cannot be ruled out with certainty until signif- 
icantly higher-resolution simulations are performed, 
it seems unlikely, because in the simulations, clouds 
with sizes down to the smallest resolved scale are often 
found even at the lowest values of pt> thus not giving 
any indication of the presence of a density-size rela- 
tion. Instead, clouds of similar average densities of- 
ten span the whole range of scales accessible to them. 
Also, the global balance between the turbulent and 
gravitational energies is mostly just a consequence 
of the presence of a self-regulated cycle of gravita- 
tional contraction, star formation, energy injection to 
the medium and dispersal, and again gravitational 
contraction, as discussed in Vazquez-Semadeni et al. 
(1995a), so that near global balance between turbu- 
lence and self-gravity is maintained at all times. 

Finally, note that here we have not discussed other 
mechanisms that have been suggested in the literature 
as responsible for the Av-R relation, such as inverse 



cascades of angular momentum, ( Henriksen & Turner 
19£4), critical thermal pressure equilibria (Chiczc 



198 7) or the contribution from the magnetic support 



(Shu 1987; Myers fc Goodman 1988; Mouschovias & 



Psoitis 1995; Gammie fc Ostriker 1995). We will ad- 



dress the role of these processes in the simulations in 
future work. 

4.3. Comparison with observations 

The lack of correlation between cloud density and 
size in the simulations appears to be in contradiction 
with the correlations found in most observational re- 



suits (Larson 1981; 


Torrellcs ct al. 1983); Dame ct al. 


1986; [Falgarone & Perault 1987 


Myers 1990|; 


WMD|). 


However, as discussed in § 


4.2.1 


we interpret the dis- 



crepancy as an effect of the limited column density 
dynamic range of the observations or of selection ef- 
fects introduced by focusing exclusively on global in- 
tensity maxima in the maps. There are some exam- 



ples of observations that have intended to avoid these 
problems ( parr 1987| ; |Lorcn 1989| ; |FPP[ ). The first two 
authors have focused on clumps within a single molec- 
ular cloud with extensive star formation, thus bet- 
ter sampling the non-gravitationally-bound turbulent 
transients. The work of FPP was specifically tailored 
towards studying low-brightness regions in molecular 
clouds. In both Loren (1989) and FPP, column densi- 
ties spanning over one and a half orders of magnitude 
are found. 

On the other hand, WMD have recently concluded 
from an analysis of 60 /im and 100 /im IRAS maps, 
that column densities of dark clouds cluster typically 
at N(H 2 ) ~ 4 xlO 21 cm~ 2 (corresponding to a typ- 
ical 100 /im opacity of r 10 o — 200 /iNepers), while 
claiming that the dynamic range of the data would 
have allowed detection of any significant variations. 
The discrepancy between these results and those of 
the present paper seemingly cannot be explained in 
terms of the limited column density dynamic range of 
their observations. 

In actuality, we believe that the results of WMD 
may be spurious and attributed to a combination of 
their selection criteria and the following effect. At 
60 /im and 100 /im, one is observing "warm" dust that 
could be coming from the "edges" of molecular clouds. 
If this "warm" dust is being heated by the visible 
photons from the galactic stellar radiation field, it is 
expected that the depth of the "warm" region will be 
of order of a few Nepers in the visible. This will lead 
to apparent constant optical depths when determined 
from observations of the "warm" dust. 

Detailed radiative transfer models have been pre- 
sented by Bernard et al. (1992), whom have already 



warned against using 100 pm surface brightness as a 
tracer of dust column density. In what follows, we 
present a simpler model that allows us to give first- 
order estimates of the scaling of the intensities in the 
IRAS 60 and 100 /im bands and the apparent dust 
opacity as a function of the true dust opacity through 
the cloud. 

Similarly to [larrett et al. (1989)| and WMD, we 
assume that the ratio of visible to 100 pm absorp- 
tion efficiency is (Qy/Qioo) = 2.0 xlO 4 . It is impor- 
tant to emphasize that the 100 /im opacity consid- 
ered in this ratio is being produced by small dust 
particles that become heated to relatively large tem- 
peratures and that are responsible for the observed 
100- and 60 /im emission. The absorption efficiency 
ratio that includes all the dust at 100 /im is about a 
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factor of 10 smaller (Chini, Krugcl & Krcysa 1986) 



Furthermore, in the range of wavelengths between 60 
and 100 fim, the absorption efficiency scales approxi- 
mately as v~ x . Using this functional dependence, we 
can write a crude approximation for the dust temper- 
ature as 

TbusT = (T c 5 r + ^ AD e^) 1 / 5 , 

where T cr and Trad are pseudotemperatures that pa- 
rameterize the heatings due to cosmic rays and radi- 
ation, respectively. (Although this approximation is 
not valid at all wavelengths, we have checked that 
using slightly different functional dependences does 
not significantly alter our results.) We assume that 
cosmic ray heating is constant for any point in the 
cloud. Specifically T cr = 10 K is the temperature 
that the dust has if heated only by cosmic rays, and 
Trad is the temperature that the dust would have 
if heated only by the stellar radiation field (at the 
edge of the cloud, where no absorption is present). 
Also, rv is the opacity in the visible, that increases 
as we get deeper into the cloud. In this simple model 
Tdust — Trad at the edge of the cloud and it tends 
to T cr for the inner regions of the cloud (fig. 12). 
The intensities at 60 /zm and 100 /im, I@o and Iioo, 
can be calculated for clouds with different values of 
Ty, as shown in fig. [l3]. This calculation is made as- 
suming that Q100/Q60 = 60/i/lOOjU. From this figure 
it is evident that for clouds exceeding a few Nepers 
in Ty, 160 and I100 do not continue growing with Ty. 



An apparent dust temperature, T, 



60/100 



can be de- 



rived from and I100 and the far-infrared opacities 
are then obtained. In fig. |lj, we plot the 100 /im 
opacity derived in this manner as a function of the 
cloud's Ty for different values of T ra d. The apparent 
100 /im opacity "saturates" at values in the order of 
70 to 100 /iNepers. Considering that along a line of 
sight one expects to intersect the front and the back 
edges of a cloud, the typical 100 /im opacities of T100 
~ 200/iNepers appear to be explained as a result of 
this effect. The few clouds with much larger values of 
T ioo (up to T100 ~ 10 4 ) reported by WMD are regions 
of strong star formation activity, which may heat the 
clouds from the inside, raising tioo- The presence of 
this effect is also consistent with the limb brightening 



observed in some clouds in the far-infrared by Sncll 



He jer fc Schloerb (1989)|and WMD, and predicted by 
the models of Bernard et al. (1992). Note that this 



effect continues to be applicable even if the clouds 
are clumpy, since it should hold at the edges of any 



density peaks, large or small, as long as they have a 
large enough column density. However, note also that 
this "saturation" effect in the determination of the 
apparent tioo will occur only if reasonably high tem- 
peratures (Tdust > 20 K) are present at the cloud's 
surface. In any case, for lower dust temperatures the 
emission at 60 and 100 /urn is very weak and unde- 
tectable in practice. One consequence of this effect is 
that in order to fully sample the dust from dark clouds 
one requires observations at longer wavelengths that 
will trace the predominant cooler dust component. 
We intend to verify this effect by comparing column 
densities obtained with different indicators in future 
work. 

The effect discussed above may clearly impose 
an upper bound to the column densities derived by 
WMD. Furthermore, concerning their selection cri- 
teria, it should be noticed that WMD define core, 
cloud and cloud complex in terms of ranges of ex- 
tinction. This directly selects against identification of 
low-column density structures. In fact, in their maps, 
small, low-extinction clouds are readily seen, but not 
classified as such precisely because of their low extinc- 
tions. Thus, the column densities of the cores studied 
by WMD are bound from above due to the "satura- 
tion" effect of the cloud edges, and from below by 
their very definition of a core, rendering their derived 
column density constancy open to question. 

5. Summary and Conclusions 

In this paper we have searched for Larson-type 
(1981) correlations and cloud mass spectrum slopes 
in the clouds generated in numerical simulations of 
the ISM, one from Paper I, a similar one with a larger 
density dynamic range, and another with higher reso- 
lution. We define a cloud as a connected set of pixels 
in the density field with values larger than an arbi- 
trary threshold pt- From the results at various dif- 
ferent times in the two simulations, we find that the 
mass spectrum has the form dN/dM oc M -1 ' 44 , 
and the velocity dispersion is related to the cloud size 
by Av oc Jl - 41 ^ 08 '^ where the errors are the standard 
deviations in the set of values including all the various 
times. The dispersion-size relation exhibits a scatter 
of about an order of magnitude, comparable to the 
scatter observed in real clouds. 

The simulated clouds do not exhibit a density-size 
relation, but instead, at all mean densities, clouds of 
sizes down to the smallest resolved scales are com- 
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monly found. This result implies that the clouds do 
not have constant column densities, but instead ex- 
hibit a range of roughly two orders of magnitude, thus 
providing strong support to the possibility that the 
observationally-derived density-size relation is an ar- 
tifact of the limited dynamic range employed by ob- 
servational surveys, and of the criteria used for se- 
lecting the clouds, as first noticed by Larson (1981), 
and then strongly argued for by Scalo (1990). Our 
results suggest that low-column density clouds do ex- 
ist in the ISM, but are systematically missed by most 
observations. In the simulations, low-column density 
clouds are turbulent transients, as in the suggestion 
by Magnani et al. (1993). Observational work that 



has used complete cloud samples (e.g., |Loren 1989[) , 
or specifically looked at structures in th e low bright- 
ness regions of molecular clouds (FTP), has indeed 
found a reported column density variability of over 
one and a half orders of magnitude, and masses well 
below the virial mass, by factors up to two orders of 
magnitude. However, note that in our simulations 
most of the mass resides in the largest clouds at each 
mean density, since there are comparable numbers of 
small and large clouds. 

The set of largest clouds at every threshold p t 
exhibits a density-size relation p ~ R a , with a = 
— 0.8±.15. These clouds appear to be close to balance 
between turbulence and self-gravity on occasions, but 
not in general. On the other hand, pure mass con- 
servation, which gives an upper limit — a < 2 for the 
two-dimensional case, together with numerical limi- 
tations of the simulations, cannot be ruled out as the 
sole causes responsible for this result. In any case, 
regardless of what the specific slope of the envelope 
turns out to be upon removal of numerical effects, 
our results suggest that the notion of a density-size 
scaling relation should probably be replaced by that 
of an "allowed" region in p-R space. Whether the 
high-column density boundary is truly a power law, 



the value of its corresponding index, and the physical 



mechanism responsible for it, arc issues that remain 
to be determined by high resolution 3D simulations. 

The result that the Av-R relation is verified in 
the simulations but the p-R relation is not, supports 
the interpretation that the former relation is a direct 
consequence of the statistical properties of the turbu- 
lence, since a turbulent energy spectrum of the form 
E(k) r~j k~ 2 , as observed in the simulations, implies 
an i? 1 / 2 scaling for the velocity dispersion. However, 
we discussed the possibility that a density-size rela- 



tion does exist, although with a scatter larger than 
the column density dynamic range spanned by the 
clouds in the simulations. Such a scatter would be 
consistent with balance between turbulence and self- 
gravity, with a scatter of one order of magnitude in 
the velocity dispersion, according to relation (^J). In 
this case, the origin of the p-R relation would remain 
unknown, and the turbulent origin of the Av-R rela- 
tion would have not as strong a support. This possi- 
bility cannot be ruled out without very high resolu- 
tion simulations in order to add at least another order 
of magnitude to the column density dynamic range. 
However, the simulations do not give any indication 
that this may be the case, since cloud sizes down to 
the smallest resolved scales are found even at the low- 
est values of pt- Besides, the column density dynamic 
range in the simulations is already larger than that of 
most observational surveys. 

We also discussed the recent results of WMD, who 
have argued in favor of a surprisingly constant column 
density (to within a factor of a few) in a large sample 
of cores studied through IRAS 60 and 100 /mi maps, 
while claiming a very large dynamic range in their ob- 
servational method. We argued that this result may 
be spurious, by presenting a model in which the warm 
"skin" of clouds and cores, which is likely to always 
have visual extinctions of order unity, puts an upper 
bound to the column densities measured by WMD. 
From below, WMD's own selection criteria eliminate 
low-column density clouds, thus mimicking a nearly 
uniform column density. 

Finally, we remark that in this paper we have lim- 
ited the discussion to the trends of the density and 
velocity dispersion with size, omitting any discussion 
of the role of the magnetic field, which is clearly im- 



portant for the dynamics of the simulations (Paper I). 
A detailed account of the energy budget in the simu- 
lations including the magnetic energy density, as well 



as surface terms is under way ( Ballesteros-Paredes & 



Vazquez-Semadeni 1996 ). 
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Now consider the turbulent crossing (or nonlinear) 
time for scale a: 



6. APPENDIX 

In this appendix we make an order of magnitude 
estimate of the effect of the numerical difussion term 
in the continuity equation. Consider the continuity 
equation, eq. (0): 



dp 
dt 



V • (pu) 



(10) 



In order to estimate the effects of difussion, we con- 
sider the extreme case when the advection term can 
be neglected, so that the mass equation becomes the 
standard difussion equation: 



- = M V P . 



(11) 



For simplicity, we consider a one-dimensional (ax- 
isymmctric) problem, whose solution is (see, e.g., 
Habbermann 1987): 



p(x,t) 



' 1 X f(x')exp(- { \ X ' r )<h'. 



Airpt 



4pt 



(12) 

where f(x') is the initial density distribution, which 
we assume to be a Gaussian, f(x') — po exp(— x' 2 /2a 2 ) 
In what follows, we will identify the width a of the 
Gaussian with the size of the cloud of interest. Equa- 
tion (O) then becomes: 



p(x,t) = pa 



t/t + 1 



1/2 / x 2 /2a 2 

e ^TJTo~Ti l (13) 



where 



_^ N ) = ^ = 9 .4lxlO- 3 f^- 
t code ) 2p \ pixels 



(14) 



is the characteristic diffusion time in units of the code, 
icode = 1-3 x 10 7 yr, and the second equality follows 
from using the value p — 0.008 (see Paper I) (note 
that the size of the integration box is 2n in the code 
units). This value of to can be made independent 
of the resolution of the simulation by choosing p such 
that /ifcjiax = est, where /c max is the maximum Fourier 
wave number in the simulation, equal to 1/2 of the 
resolution. 



t nl ~ a/ui, 

where it; is the turbulent speed associated to the scale 
I. Note that we are allowing for the possibility that 
clouds of size a are generated by turbule nt s treams 
of different, typically larger, size I (see §hl). This 



velocity can be estimated from the turbulent energy 
spectrum Elk) as: 



1 2 
—uf = 

2 1 



Elk) dk, 



2-rr/l 



that is, ui is the root mean square energy per unit 
mass in scales smaller than I. Using E(k) = O.lfc -2 , 
as indicated by fig. [HI one obtains: 



14.05 



Z/pixels \ 1 / 2 
i rcs /pixels/ 



(15) 



where n rcs is the number of pixels (i.e., the resolution) 
per spatial dimension of the simulation. From equa- 
tions (Q) and (|l5|), we can now compare the diffusion 
and nonlinear times, in order to determine the scale 
at which they are equal. Assuming that I = ma, we 
obtain: 

/ l eq \ 16.33 
\ pixels J m 1 / 2 ' 

and thus, clouds with sizes lower than ~ 16 pixels are 
dominated by the numerical mass-difussion. There- 
fore, it is important to assess the effect of mass dif- 
fusion in the results of §||. To do this, we take the 
extreme position that for time durations At < tnl, 
a Gaussian cloud is affected exclusively by diffusion, 
and compute the fractional variation of its central 
density over tnl ■ In fig- [l5| we show the evolution of 
p(0,t)/po given by equation dl3| ) over tnl for a cloud 
of size a — 4 pixels. Note that this is a worst-case es- 
timate, since on the one hand, in the plots below we 
have discarded clouds with sizes < 4 pixels in order 



to avoid viscosity effects (see 64.1), which is most af- 
fected by diffusion, while, on the other hand, we have 
taken m = 1. 

As we can see from fig. the density decreases 
by a factor <J 5. Typically, thus, we can expect our 
density data to change by factors smaller than half an 
order of magnitude. 

Using the above results, we can produce a "cor- 
rected" density-size plot with the estimated "true" 
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densities p C orr calculated as: 

Peon ~ Pdata(l + ^NL Ad) ^ ■ 



(17) 



where pdata are the raw cloud mean densities as pro- 
duced by the simulations, and the "true" sizes are 
obtained assuming mass conservation: 



Rcovr 



Rda 



(I + tnlAo) 1 ^' 



(18) 



The "corrected" density-size plots are shown in figs. 
p"6| and [I?] for Run 28bis and the high-resolution Run 
28.800, respectively. We see that, although slightly 
pushed closer to a correlation in the case of Run 28 
bis, the general trend of these plots still supports the 
main conclusion from §||, namely that clouds of sizes 
down to the resolution are seen at all mean densities. 
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Fig. 1. — Density (contours) and velocity (arrows) 
fields for Run 28 at t = 6.6 x 10 7 yr. The contours 
are spaced logarithmically, with increments of 0.3125 
in logp. The minimum density in this plot is 0.04 
cm~ 3 , and the maximum density is 40 cm~ 3 . 




Fig. 3. — Logarithmic mass distributions for a) Run 
28 at t = 6.6 x 10 7 yr and b) Run 28bis at t = 7.15 X 
10 7 yr. The lines show least-squares fits to data bins 
with logM > 0.8, and have slopes —0.51 in a) and 
-0.43 in b). 




Fig. 2. — Regions ("clouds") with densities larger 
than an arbitrary threshold pt in the density field of 
Run 28 at t = 6.6 x 10 7 yr. a) p t = 4. b) p t = 16. 
Note that small clouds with high densities are nested 
within larger complexes, but small clouds with low 
densities can be seen in a) as well. 





Fig. 4. — Velocity dispersion vs. size for all clouds in 
a) Run 28 at t = 6.6 x 10 7 yr and b) Run 28bis at 
t = 7.15 x 10 7 yr. Clouds obtained with each value 
of pt are shown with a specific symbol as indicated. 
The solid lines show least-squares fits, with the slopes 
and correlation coefficients indicated. For reference, a 
1/2 slope is indicated by the dotted lines. The cloud 
labels in a) and b) are respectively the same as in a) 
and b) of figs. |k and b. 
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Fig. 5. — Average cloud density us. size for for all 
clouds in a) Run 28 at t = 6.6 x 10 7 yr and b) Run 
28bis at t = 7.15 x 10 7 yr. The lines show a p oc R^ 1 
power law. In both cases, the largest cloud at each 
Pt lies near the p oc Rr 1 line, although other times 
tend to show shallower envelopes (see fig. ^). The full 
ensemble of clouds does not show any trend of (p) with 
R. The cloud labels in a) and b) are respectively the 
same as in a) and b) of figs. ||a and b. 



Fig. 6. — Average density vs. size for the largest 
cloud at each value of pt at six different times in Run 
28 (bottom curves) and three times in Run 18bis (top 
curves, labeled "b"). Each point denotes the largest 
cloud at one value of p t . Different values of pt for the 
same time in the simulation are joined by dotted lines 
and indicated with the same symbol. The various 
curves are displaced by increments of 0.2 in logp for 
clarity. The numbers next to each curve show the time 
into the evolution of the run to which it corresponds, 
in units of 10 7 yr. Note that time t = 6.6 x 10 7 yr for 
Run 28 and time t = 7.15 x 10 7 yr for Run 28bis have 
some of the steepest slopes. The straight line shows 
a p oc power law. The average slope for all times 
is —0.81, with a standard deviation of 0.16. Because 
of the "compression" at high densities due to various 
numerical effects, only clouds with logi? > 1 (Run 
28bis) or logi? > 1.5 (Run 28bis) were considered in 
the calculation of the average slope (see text). 
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Fig. 7. — Logarithmic grayscale image of the density 
field of run 28.800 at time t = 7.15 X 10 7 yr. This run 
is similar to run 28bis, except that it has larger reso- 
lution (800 x 800), and that star formation is turned 
off at time t — 6.9 x 10 7 yr. The density extrema are 
Pmax = 129, and p min = 3.4 x 10~ 2 . 
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Fig. 9. — "Familiy tree" of a few selected clouds from 
fig. ||b. The largest cloud is seen to branch off to 
both clouds defining a Larson-type relation as well 
as to clouds away from it (dotted lines). A similar 
branching pattern is observed for clouds of lower col- 
umn densities (dashed lines). 



Fig. 8. — Mean density-size plot, equivalent to figs. 
||a and b, but for the high-resolution Run 28.800. A 
large range of sizes, in this case spanning over two 
and a half orders of magnitude, is again seen at the 
lowest values of the average cloud density. 
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Fig. 10. — Turbulent energy spectra for Run 28 at 
t = 6.6 x 10 7 yr. Solid line: incompressible com- 
ponent of the spectrum. Dashed line: compressible 
component. The straight line shows a k~ 2 power law, 
characteristic of shocks, which implies an Rr 1 / 2 de- 
pendence of the velocity dispersion. 
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DUST OPACITY IN THE VISIBLE (NEPERS) 

Fig. 12. - Dust temperature, Tdust, as a function 
of dust opacity in the visible for Trad = 25 K. Near 
the edge of the cloud (left side) Tdust tends to Trad , 
while toward the inner regions of the cloud (right side) 
heating by radiation becomes negligible and Tdust 
tends to the value provided by cosmic ray heating (10 
K in our model). 



Fig. 11. — Evolution of the total turbulent kinetic 
(solid line) and the gravitational (dotted line) energies 
per unit mass in code units for run 28 over the last 
half of its evolution. The two energies are very close 
to equipartition at all times. 
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Fig. 13. — Intensity of cloud emission at 100 /xm and 
60 fim as a function of dust opacity in the visible, for 
Trad=25 K. Note that the growth of the intensity 
"saturates" above a few Nepers of dust opacity in the 
visible. The reason for this effect is that beyond a 
few Nepers of dust opacity in the visible there is no 
significant radiation heating and the dust becomes too 
cold to emit significantly at 100 fim and 60 fim. 



Fig. 14. — Apparent warm dust opacity at 100 /zm 
(derived from the intensities at 100 /mi and 60 fim), 
as a function of dust opacity in the visible for different 
values of Trad- Note the "constancy" in the appar- 
ent dust opacity at 100 fim once the cloud exceeds a 
few Nepers in dust opacity in the visible, practically 
independently of Trad • 
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Fig. 15. — Fractional change of the central density 
p(0,t)/po of a cloud of size 4 pixels. The fractional 
change shown here is an upper bound to the correc- 
tions applied to the simulated data, since we have 
considered a cloud with size smaller than the mini- 
mum cloud size retained, and to = 1 in eq. (|l6|). 
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Fig. 17.— Same as figure ||, but for a Run 28.800. 
Here, the sizes of the smallest-density clouds vary by 
nearly 2.5 orders of magnitude. Note that, although 
only clouds with sizes larger than 4 pixels are re- 
tained, the smallest clouds in this figure have sizes 
smaller than that because of the size "correction", 
eq. ©. 



Fig. 16. — "Corrected" density-size plot for Run 28bis 
using the estimates for the 'true" density and size of 
the clouds given by equations (|l7|) and ([[si). Note 
that a range of roughly two orders of magnitude in 
size remains for the lowest mean density clouds. 
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